Manipulating energy and spin currents in nonequilibrium systems of interacting qubits 
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We consider generic interacting chain of qubits, which are coupled at the edges to baths of fixed 
polarizations. We can determine the nonequilibrium steady states, described by the fixed point 
of the Lindblad Master Equation. Under rather general assumptions about local pumping and 
interactions, symmetries of the reduced density matrix are revealed. The symmetries drastically 
restrict the form of the steady density matrices in such a way that an exponentially large subset of 
one-point and many-point correlation functions are found to vanish. As an example we show how 
in a Heisenberg spin chain a suitable choice of the baths can completely switch off either the spin 
or the energy current, or both of them, despite the presence of large boundary gradients. 

PACS numbers: 75.10.Pq, 03.65.Yz, 05.60.Gg, 05.70.Ln 



I. INTRODUCTION 

The impressive progress in experimental manipulation 
of nanowires and quantum dots makes it possible to in- 
vestigate quantum systems consisting of a few quantum 
dots or quantum bits. On the other hand, manipula- 
tions/operations on a single quantum bit (e.g., applica- 
tion of a quantum gate) are at the basis of the func- 
tioning of any elementary block of a quantum computing 
device. However, a theoretical understanding of micro- 
scopic quantum systems out of equilibrium (e.g. under 
constant pumping or continuous measurement by a quan- 
tum probe) is far from being complete, apart from simple 
cases like a single two-level system or a quantum har- 
monic oscillator under external pumping or in contact 
with reservoirs 0, Q ■ On the other hand, using the dis- 
sipative dynamics for the preparation of quantum states 
with required properties is becoming a promising field of 
research [1, HI]. In this respect, the role of sufficiently 
simple spatially extended systems, amenable to both an- 
alytic and numerical investigations, becomes important. 

Unlike a quantum system evolving coherently, whose 
evolution at any time depends on the initial state, a 
quantum system with pumping tends to a steady state, 
independently of the initial conditions. This allows to 
manipulate the steady state of the system by a suitable 
choice of the pumping. In this paper we first illustrate 
how to accomplish such a task for a general system of 
interacting quantum spins. In particular we consider 
symmetries of the density matrix that impose rigid con- 
straints on the properties of the nonequilibrium steady 
state, thus entailing, for instance, exact vanishing of cu- 
mulative correlation functions, like certain components of 
the structure factor. Then, we specialize our analysis to 
the one dimensional driven XXZ chain of quantum spins 
in the presence of pumping applied at the edges. We show 
that, for particular realizations of the pumping and irre- 
spectively of the system size, one can switch on and off 
the spin and/or the heat currents in the nonequilibrium 



steady-state. Such an approach unveils new interesting 
perspectives in the study of driven spin chain models. 
They have been mainly investigated to understand un- 
der which conditions the spin and the energy currents in 
the steady state exhibit anomalous or ballistic behavior 
In this perspective, we want to point out that the 
results presented in this paper hold independently of the 
anomalous or ballistic features of transport [8| and even 
in the absence of integrability. 

In Section [Tf] we present the general model of Lindblad 
master equation and discuss how symmetries may affect 
its properties. Specific forms of Lindblad operators act- 
ing on quantum spin models are introduced in Section Hill 
In particular, we devote special attention to parity selec- 
tion rules in Section HVl Section [V] deals with the special 
case of the Lindblad dynamics for the one dimensional 
XXZ model of qubits. It has been chosen to illustrate 
in detail how our findings apply to a simple model, that 
has been the object of intense resent research. Specifi- 
cally, in subsection lV Al we discuss the kind of symmetries 
emerging in this model for Lindblad operators of target 
type, while subsection IV Bl is devoted to the analysis of 
spin and energy conductance in the presence of gradients. 
Conclusions and perspectives are contained in Scction fVTl 



II. 



LINDBLAD MASTER EQUATION AND ITS 
SYMMETRIES 



We consider the quantum Master equation in the Lind- 
blad form 0, @ , 

(i) 

where we have set h = 1; p is the reduced density matrix, 
H is the Hamiltonian of the system and Lm is the Lind- 
blad operator. It is easy to verify that JjTr(p) = 0, p + = 
p, ^Tr{p 2 ) 7^ 0. The first two relations are necessary for 
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interpreting p as a density matrix, with Tr{p) = 1, while 
the latter implies that we deal with an open system: an 
initially pure state p = \4*)(4>\ will not remain pure in the 
course of time. The Lindblad equation is the most gen- 
eral equation of motion for the reduced density matrix, 
which assures positivity and conserves the trace 

The non-unitary part of the Lindblad Master equation 
(LME) makes the dynamics irreversible. In the course 
of time, any initial state will relax to a steady state, de- 
scribed by the time-independent solution of the Master 
equation. Our purpose is to reveal generic properties of 
this stationary state, namely constraints imposed by the 
symmetries. In the simplest case, such constraints take 
the form of selection rules, according to which a subset 
of the matrix elements of the reduced density matrix has 
to vanish. 

Let us denote the right-hand side of LME as C[p]. 
If it is invariant under a unitary transformation U , i.e. 
C[UpU + ] = UC[p]U+, it follows that p(t) = Up{t)U+ is 
a new solution of LME. In general the solutions p{t) and 
pit) describe different time evolutions. If the nonequi- 
librium steady state (NESS) of the system pness — 
lim^oo p{t) is unique [10(, then the trajectories p(t) and 
p(t) eventually converge in time to the same asymptotic 
solution, limt^-oo p{t) = limt^-oo p(t) = Pness- Accord- 
ingly, the unique steady state of the system has to be 
invariant under the transformation U, 

Pness = Up N ESsU + . (2) 

This implies also that the expectation value of any phys- 
ical observable of interest /, measured in the steady 
state, (/) = Tr (^fpNESs), has to satisfy the relations 

(/) = Tr (fUp N ESsU+) = Tr (u+fUp N ESs) ■ 
In particular, if / changes sign under the action of U, 
U + fU = -f, (3) 

it follows that 

(/> = Tr (U+fU pness) 

= -Tr (/ pness) = -(f) -> (/) = . 

On the contrary, if / is invariant under the action of U, 
i.e. 

U+fU = f , 
no consequences for (/) can be drawn, 

(/> = Tr \ U + fU pness) = Tr (fp N Ess 



III. 



LINDBLAD OPERATORS 



For instance, let us suppose that a two-level open 
quantum system, described by a generic density matrix 
p = |/+ \ ^Z a {< ja )< Ta {& a being Pauli matrices), is in- 
variant under the transformation U = a z , and has a 
unique steady state. Since a z a x a z = —a x , a z a y a z — 
—<r v and <7 z a z a z — a z , from ^ we can conclude that 
(<j x ) = (<r y ) = 0, while (a z ) is invariant under the action 
off/. 



Here we focus on the most commonly used types of 
Lindblad operators: (i) those targeting a given value of 
observables (tTp); (ii) those responsible for dephasing, or 
decoherence, whose action favors the evolution to a clas- 
sical state; (iii) the non-local variants of the previous 
cases. 

Lindblad operators targeting a given value of z— spin 
projection afarget a ^ a chosen site p have the form of 
creation-annihilation operators, 



L[ p) = a{a x p + i<r») = aa+ 
Lf = (3(a x p - »o*) = fia- 



(4) 
(5) 



The equations of motion for the expectation values of 
the operators Cp , c p , a z read 



d(<J z p ) 

dt 
d(<r x ) 

dt 

dt 



H(a z p )-T z ((a; 
U{a x p )-T x (a x p 



t) 



(6) 



where we introduce the shorthand notation H{f p ) = 
-iTr(f p [H,p\), and 



T z =4(a 2 +(3 2 ); a z algct 



/3 2 



(7) 
(8) 



If the coupling constant T z is sufficiently large with 
respect to the norm of the Hamiltonian in equations 
([6|), then % can be neglected and the averages {u p (t)} 
converge, after some relaxation time of order 1/r^, to 
their "targeted" values crf arKOt = 0,of arget = 0, of argct = 



■0 1 



target 

as follows 



K(o) 



u target / e 



(9) 



Lindblad operators of type §5§ naturally appear in the 
problem of an atom interacting with a quantized radia- 
tion field [11] , in spin chains coupled to a bath of fixed 
polarization, in electron paramagnetic resonance experi- 
ments and in studies of quantum transport [l2| - [2l| . 

Lindblad operators targeting a given value of x— spin 
projection or y— spin projection are given by the appro- 
priate cyclic rotation of the operators (O, i-e. 

Vi = a(a v p +ia z p ), V 2 = flip* - ia z p ) (10) 

with target crf argct = (a 2 — /3 2 ) / (a 2 + /3 2 ) at site p and 

W 1 =u{a%+%(%), W 2 = v{a z p -ia x p ) (11) 

with target cr^ argct — (u 2 — v 2 ) / (w 2 + v 2 ) at site p. The 
relaxation to the targeted values can be described in com- 
plete analogy with (U])-©, see (22J. 
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Dephasing Lindblad operators have the form 



L 



(p) 

deph 



(12) 



and describe the presence of a dephasing noise in the 
dynamics, which causes decoherence (i.e., vanishing of 
non-diagonal density matrix elements) with rate 7 P at 
site p. The local equations of motion, analogous to (JB]), 
have the simple form 



dt 



dt 



= 0. 



(13) 



The coupling to a heat bath is often modelled by the 
application of the dephasing Lindblad operators (|T2|) at 
all sites of the system. 

Lindblad operators, acting on more than one site, in- 
troduce incoherent non-local processes. For instance the 
Lindblad operator 



L 



(p. g) 

hopp 



VlVpVq 



(14) 



describes incoherent spin flips between sites p and q [2c 



IV. PARITY SYMMETRY SELECTION RULES 
FOR LINDBLAD DYNAMICS 

We consider a generic open system of qubits, with in- 
ternal pair interactions, described by the Lindblad Mas- 
ter equation (pj, with the Hamiltonian 



N 



H= J2 Jx{k,m)°l° X m + Mk,m)atal (15) 

h,m— 1 

N 

+ J z (k, m)a z k <7 z m + h k°l , 
fe=i 

where J a (k,m) are couplings between qubits k, m, and 
hk are local magnetic fields. We do not impose any re- 
striction on the space dimension or on the geometry, but 
we just assume the connectivity of the graph. The sys- 
tem can be subject to external pumping and/or external 
noise, modelled by the z-polarization targeting operators 
([5]), dephasing Lindblad operators (fl~2|) and incoherent 
hoppings (|14[) , with all of these operators acting on an ar- 
bitrary subset of sites. Then, if the steady state is unique, 
the transformation fi z = {(j z )® N = a z <g> a z <g> ... ® a z , 
fij = Q z identifies a symmetry of the Master equation, 
thus yielding the relation 



PNESS = ^zPNESS^z 



(16) 



Indeed, Hamiltonian (fl5|) as well as dephasing Lindblad 
operators (fl"2"|) and incoherent hoppings (fl"4")l are invari- 
ant under f2 z , while z— polarization targeting operators 
§5§ change sign under its action. Since the Lindblad part 
of the evolution equation is quadratic in L m , LME is in- 
variant under fl z . The symmetry (|16|) is known in spin 



models as parity symmetry (P-symmetry). In the special 
case where the total set of Lindblad operators contains 
either only dephasing Lindblad operators (fT2|) or only z- 
polarization targeting Lindblad operators ([5]), LME has 
one further symmetry, the PT- symmetry, which has in- 
teresting consequences on the full spectrum of the Lind- 
blad superoperator [24] . In our general setting, where 
the dephasing and the polarization targeting Lindblad 
operators are mixed, the quantum Lindblad dynamics is 
not PT- invariant. Another remark concerns our crucial 
assumption of a uniqueness of the steady state. The ex- 
istence and uniqueness of the steady state is guaranteed 
by the completeness of the algebra, generated by the set 
of operators {H, L m , L+ } under multiplication and addi- 
tion |10j | , and it is verified straightforwardly as in [25j , for 
any choice of a set of the Lindblad operators {L m }, pro- 
vided the set contains at least one polarization targeting 
operator ^ or ([5]). 

The P-symmetry yields severe limitations on the 
nonequilibrium steady state pness- We indicate with 
P l jih " % jN ^ ne matrix element of pness m the natural 
basis, labelled by indexes i\,i2,--3N which take values 
— 1,1. Let us calculate how the matrix element p£*? "[ *^ 
changes under the action of the P-symmetry operator 
£l z . One obtains 



((a z r»p(a z ) 



.IN 
■3N 



N 

n 

rn — 1 



The factor K 



N 

n 



ij m may only take values 1 and —1. 



If K = 1, the P-symmetry does not yield any constraint 
on p. Conversely, if K = —1, from (|L6|) it follows that 
the corresponding matrix element vanishes, 



iV 



3132 



■\ N =0 

■■3n 



* n 



(17) 



m — 1 



We call this condition a parity selection rule (PSR). If 
(JTTJ) holds, simple analysis shows that each row and each 
column of p contains 2^^ zero entries. For instance, 
the first row of the density matrix for 3 sites /0„^L con- 
tains four null elements: p}} 1 !! = p_\ \ i = p\-li = 
p l _\ lj _ 1 = 0. More generally, it can be easily real- 
ized that by a suitable reshuffling of rows/columns the 



2 density matrix satisfying p7|) can be brought 
into a block-diagonal form, with two blocks of equal size 

2 N-1 x 2 JV-1_ 

An important feature of PSR is that any subsystem 
made of n qubits, described by the reduced density ma- 
trix p( n ) — TrN-nP, keeps the same symmetry, as it can 
be easily verified, 



Pi V 1 ' 2 ' 



for J] 



For n — 2 the explicit form of the generic density matrix 
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(a b' 

c d 

d* ci 

\6* aiy 



(18) 



States of the form (118)) are well known in information 
theory as A-states and they are subject of intensive in- 
vestigation, (see e.g. HH). 

Finally, we want to point out that, in terms of physical 
observables, PSR ([T7l) entails vanishing of a set of exper- 
imentally measurable quantities, like many-point corre- 
ct), and structure factors 

ik(m—n) I „ol „P 



lation functions (ct^ct^ 



(«+i>, namely 







(19) 



discussed in [24J , in this case there exist two further sym- 
metries 



PNESS — ^xRPNESsR^x 



Pness — £l y Rp 



NESS 



(25) 



(26) 



= (C(g>....(g)B<S)A)R is a left-right 
® N . We shall see that the above 



where R(AigiB®...igiC) 
reflection and fl a = (<r a 
symmetries impose restrictions on transport properties 
in the driven XXZ chain. The transport properties are 
governed by the spin and the energy current operators, 
jn i7n and J„ , which are defined by the lattice continuity 

equations ^ t a z n = J n -i, n -3n,n+i, j- t h n .n+\ = J„ 
where 



fE 



Jn,m — z '\ u n u m 



T v a x ) 



(27) 



(20) 



and 



+ A(j„_i,„CT^ +1 +0-n-lin,n+l)- (28) 



( CT n CT mi <7 ro 2 "- <7 mJ = (^XA'^mJ = ( 21 ) 



S xz (k) = S yz (k) = S z:c (fc) = S*»(fc) = o. 



(22) 



(23) 



In Refs [27j and [HI the structure factors S al3 (k, A) were 
proposed as entanglement witnesses: they are measurable 
quantities in neutron scattering experiments. 



V. ID DRIVEN XXZ MODEL 
A. Parity symmetry and energy current 

The general properties described in the previous sec- 
tion can be specialized to the study of further symmetries 
emerging in one-dimensional driven spin chain models 
with pumping acting at the edges. In these cases Lind- 
blad operators create effective boundary gradients. A 
commonly studied setup (see [l2Tj20| ) is the XXZ spin 
chain, whose Hamiltonian reads 



H = 



N-l 

E 

k=l 



a x k a x k+1 +alat +1 +Aa z k a z k+1 . (24) 



There are four Lindblad operators, £1.2 = yT(l 
and L 3)4 = v/T(l 



P)vf 



fi)af 
1 < P < 1 



> N , parametrized by 
and F > 0, that target at the chain edges spin configura- 
tions with equal amplitude and opposite sign, (af) [i. 
((T N ) —> —fi, see Eq. ([7]). In addition to the P-symmetry, 
discussed in the previous section and the PT-symmetry 



It can be easily checked that the energy current operator 
J,f changes sign under the action of any of the sym- 
metries (ESI),®, n x Rj^m x = n v Rj^m y = -j^, 

thus implying that in the steady state ( J r f ) = for any 
system size. On the other hand, the magnetization cur- 
rent (|27p is invariant under the above transformation, 
Q x Rj n Rtt x = Q, y Rj n RQ y — j n and therefore it is al- 
lowed to flow. So magnetization current can flow and 
the energy current is suppressed completely, despite the 
presence of boundary gradients. 

As a validation of our results, we have integrated nu- 
merically equation (JlJ that contains, in addition to the 
operators L\- L4, also the operator V — p{a v N + «cr^), 
acting at site i = N. The presence of this additional op- 
erator V in (QJ invalidates the P-symmetry (IT6|) . In Fig. 
[1] we plot various one- and two-point correlations as a 
function of the amplitude p of the P-symmetry breaking 
operator V: as expected they are found to vanish only 
for p = 0. 

A diffe rent choi ce of the Lindblad ope rators, n amely 
Li, 2 = v / T(iTm)(^ ±io\) and L 3A = ^T(l - f i)(a z N ± 
i(rff) amounts to set a boundary twisting gradient in the 
Ay-plane: the P-symmetry (fi"6|) is violated, but other 
symmetries appear, predicting a phenomenon of a sign 
alternation of the magnetization current with the system 
size (see Q, [29j|). 



B. Spin and thermal conductance with boundary 
gradients 

There is a large interest in studying the conductance 
in low-dimensional materials, due to the rich and often 
counterintuitive features they exhibit. In this section we 
discuss in full generality how one can switch on and off 
the magnetization and the energy currents by a suitable 
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Figure 1: Some observables characterizing the steady state 
solution of the Lindblad Master equation {1} with an addi- 
tional Lindblad operator V = p{cr v N + io z N ) versus the ampli- 
tude p. Thin, broken and thick lines correspond to one-point 
correlations (a„) , (o - ^) and spin current respectively. 
For p — the symmetry 1)170 is valid and all observables 
vanish. The adopted parameter values are N = 4, Jx 
/v = 1, J z = —0.6. The Lindblad operators are L\ = 2af , 
L 2 = L 3 = 0,1/4 = \/2o-^, V = p(ct^ + io-^r). 



choice of boundary reservoirs , i.e. Lindblad operators, 
acting on the XXZ spin chain ([24]) . 

Let us couple the XXZ chain at the boundaries to 
baths of constant (but different) magnetizations, so that 
the time evolution of the state becomes dissipative and 
is described by LME 



dp 
Of 



-£ = -i[H,p]+T(£ L [p} + C R [p]), 



(29) 



where H is the XXZ Hamiltonian of the open XXZ 
chain with anisotropy A (see Eq. (pM|) ) . £l[p] and 
are Lindblad dissipators acting on the leftmost 
(k = 1) and on the rightmost (fc ~ N) boundary spins, 
while r denotes the interaction rate with the dissipators. 
By choosing different parameter values of the bound- 
ary Lindblad dissipators C l [p] and Cr [p] , spin gradients 
can be introduced. In particular, the Lindblad dissipa- 
tors target spin polarizations at site 1 and at site N, 
described by the one-site density matrices pl and pn 
satisfying Cl[pl] — 0, Cr{pr] — 0, respectively. For 
sufficiently large values of T the reduced density matrix 
of the system p(t) evolves in time towards a nonequi- 
librium steady state density matrix, pnesSj such that 
Tr 2 ,3,...NPNESS -s- PL and Tri, 2 ,...N-\PNESS PR- 
Let us choose the following Lindblad dissipators: for 



S X L,R,S y L ,R 
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Figure 2: Targeted spin projections at the left and at the right 
boundary s x R (full line), s L (full bold line) and s v R (dashed 
line), s\ (dashed bold line) versus a, from (|32)l and (|33l) . with 



the left boundary, 

C L [p] = AV 2 pV+-- {p, V 2 + V 2 }+aW lP W+-- {p, W+W,} 

(30) 

where V 2 = o\ — i<rf, W\ = of + icrf , and for the right 
boundary, 

Cu[p] = V zP V+- 1 - {p^V+V^+AaW+pWi-^- {p,W+W A } 

(31) 

where V 3 = a v N + ia z Nl W^ = a z N — ia x N . We assume also 
the following conditions, 4 ^ 1 and < a < 1. 

It is straightforward to check that £l[pl] = where 

pl = f + \s x L a x + \s\oy + \s z L ay with 



2A 



2a 



2A + a 



> °L 



A + 2a 



; si = (32) 



The set {s L , s v L ,s L } are targeted spin components at the 
left boundary. At the right boundary we have to choose 

2 „ 2aA 



2 + aA 



l + 2aA' 



4 = (33) 



in order to obtain £r[pr] — 0. The graphs of the tar- 
geted spin components at the left and at the right bound- 
aries for A = 2 are shown in Fig(5] The left bound- 
ary values arc different from the right boundary values, 
(s L (a),s v L (a)) 7^ (s^,(a), s v R (a)), entailing a presence of 
boundary gradients for any value of a. 

Due to the Heisenberg exchange interaction, one might 
expect that the presence of a spin gradient yields non- 
vanishing spin and heat currents, given by the Fourier 
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law 



}■■' 



As^ 
Al 
As? 
Al 



where As^/Al 



(4 



s l)/( n - !)■ As shown in FigH 



for all values of < a < 1, targeted spin components at 
the left and right boundary are different, thus yielding 
an applied gradient. However, we find that for a = 
the steady spin current is identically zero, (j) = , while 
(J E ) ^ 0. On the other hand, for a = 1 we obtain the 
opposite scenario, i.e. (J E ) = and (j) ^ 0. In fact, for 
a = the stationary solution of the Lindblad equation 
Pness is invariant under the following transformation, 



Pness = QxPNESS^x , 

where fl x — (a x )® N . Analogously, for a 
invariant under the transformation 



(34) 

1) PNESS IS 



PNESS = QxUrotRpNESsRU rot Q x 



(35) 



where R is again the left- right reflection, (A ® B eg) ... ® 
C) = {C®....®B® A)R, and the diagonal matrix U rot = 
diag(l,i)® N is a rotation in XY plane: U ro to' 3 ^lJ^ ot — 



r« IT t rryU + 

'hi Ureter. U rr .+ 



-a*- The Hamiltonian part of LME, 



— i[H, p], is also invariant under both transformations, 
while for the Lindblad part the symmetries are satisfied 
due to the specific forms of C E [p] and Cr [p] for a = 
and a = 1. 

Case a = 0. Making use of the symmetry (13"4")) and 
of the properties of the Pauli matrices, we obtain the 
following expressions for the magnetization and for the 
energy currents (in what follows we use the shorthand 
notations j and J E for these quantities), 

j = Tr(p NESS j) = -Tr(E x p NE ss^xj) = -j 



J E = Tr(p NESS J E ) = Tr(Z x p NESS X x j E ) = J E 

The first one of these relations implies j = 0, while no 
restrictions are imposed for J E . 

Case a = 1. We find the opposite situation: the en- 
ergy current J E under the transformation ([35]) changes 
sign, while no restrictions are imposed for the magneti- 
zation current j. We conclude that J E = 0. 

Case < a < 1. For any intermediate value of a, 
neither (|34|) nor ([35]) are satisfied. Consequently, both 
magnetization and energy currents are allowed to flow. 

In order to check our findings, we have obtained nu- 
merical solutions of LME for small sizes N and dif- 
ferent values of the parameter A ^ 0. In all cases we 
find complete agreement with theoretical predictions. A 
typical case is illustrated in Fig|3j For A = 1, in addition 
we find that both J E =Q A=1 = and j a =o,A=i — 0, while 
only j = is predicted by the symmetry (|34[) . Looking 




Figure 3: Steady magnetization and energy currents as a 
function of a, from the numerical solution of the Lindblad 
equation (|29[1 . The adopted parameter values are A = 2, N = 
5, A = 1,T = 0.5. 



for an explanation, we readily find another symmetry of 
valid for A = 1 and a = 0: 



Pness — QxRpnessR^x 

which explains why also J E =0 a—i — 0- I n fact, it 
can be easily checked that under this symmetry the en- 
ergy current operator changes sign, Tr(pMESsJ E ) — 
-Tr(Z x R PNES sRZ x J E ). 



VI. CONCLUSIONS 

Transport properties of quantum systems can exhibit 
unexpected features if the noncquilibrium steady state 
has to obey certain symmetry properties. Various exam- 
ples have been discussed in this paper for models of inter- 
acting systems of qubits, subject to local pumping mech- 
anisms from specific Lindblad operators. We have first 
introduced different classes of these operators as spin- 
targeting and dephasing ones. Then, we have discussed 
the kind of symmetries they impose to the Lindblad mas- 
ter equation and to the corresponding noncquilibrium 
steady state. The important role played by parity sym- 
metry selection rules has been illustrated for a general 
Hamiltonian model. These considerations have been also 
specialized to the XXZ spin chain model. We have shown 
that spin and energy currents can be suitably regulated 
by acting on the symmetries of the NESS through the 
parameters of the Lindblad operators. In particular, we 
find that both currents can vanish, even in the presence 
of finite applied gradients. 

We have to point out that all the results reported in 
this manuscript rely on the basic assumption of unique- 
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ness of the steady state solution of the Lindblad Master 
equation. Such a property applies to all the examples 
considered in this paper. An explicit check of this prop- 
erty can be performed by using the completeness criterion 
of the algebra generated by the Hamiltonian and by the 
Lindblad operators [Iol |. Once the uniqueness is estab- 
lished, the nonequilibrium steady state is invariant under 
all the symmetries of the Lindblad master equation. In 
fact, any violation of a symmetry yields the existence of 
at least a one-parameter family of steady state solutions 
as a direct consequence of the linearity of the equation 



In a general perspective we can affirm that the proper- 
ties of the steady states analyzed in this manuscript can 
be viewed as a first achievement in the exploration of new 
interesting features of the quantum Master equation. In 
Sec|V]we have also shown an explicit example of how the 
vanishing of a current signals the presence of additional 
symmetries. 
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